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Abstract 

The influence of gluon and Goldstone boson induced tensor interactions on the dibaryon masses 
and D-wave decay widths has been studied in the quark delocahzation, color screening model. 
The effective S-D wave transition interactions induced by gluon and Goldstone boson exchanges 
decrease rapidly with increasing strangeness of the channel. The tensor contribution of K and t] 
mesons is negligible in this model. There is no six-quark state in the light flavor world studied 
so far that can become bound by means of these tensor interactions besides the deuteron. The 
partial D-wave decay widths of the IJ^ = ^'^'^ Nil state to spin and 1 AH flnal states are 12.0 
keV and 21.9 keV respectively. This is a very narrow dibaryon resonance that might be detectable 
in relativistic heavy ion reactions by existing RHIC detectors through the reconstruction of the 
vertex mass of the decay product AH and by the COMPAS detector at CERN or at JHF in Japan 
and the FAIR project in Germany in the future. 

PACS numbers: 12.39.-x, 14.20.Pt, 13.75.Cs 



I. INTRODUCTION 

There might be two kinds of dibaryon|l|, |2|. One is the loosely bound type consisting 
of two octet baryons; the deuteron is a typical example. The others are tightly bound; the 
H particle had been predicted to be such a six quark state although later calculations cast 
doubt on itJ3l,l4|. Instead, a non-strange IJ^ = 03^ d* and a strangeness -6 IJ^ = 00^ di-Jl 
have been predicted to be tightly bound six quark states, which are formed from decuplet 
baryons 2, S la III la Isl • The strangeness -3 IJ'^ = \2'^ Ni7 has also been predicted to be of 
the tightly bound type^ lid, |ll|. 

The tensor interaction due to vr exchange plays a vital role in the formation of loosely 
bound deuteron. In the d* case the tensor interaction contribution to its mass is minor but 



is critical for its D-wave decay to the NN final statelo. There are other near threshold and 

nn 

deeply bound dibaryon candidates found in two systematic quark model calculations!^, |8||. 
This naturally raises the question as to whether or not the tensor interaction adds sufficient 
strength to bind these other near threshold states to become strong interaction stable, as in 
the deuteron case? Conversely, is the tensor interaction weak enough to leave the hig h spin, 
deeply bound states as narrow dibaryon resonances, as was shown in the d* case(l2|? 

The present calculation is aimed at answering these two questions for the dibaryon candi- 
dates in the u,d,s three flavor world within the extended quark delocalization, color screening 
model (QDCSM). Our results show that both the effective S-D wave transition interactions 
due to gluon and vr exchanges decrease rapidly with increasing strangeness, and that the 
tensor contributions of K and rj mesons are negligible after a short range truncation. Al- 
together, the tensor contributions are not strong enough to bind other near threshold six 
quark states, such as the SIJ^ = — 401+HS, to become strong interaction stable with the 
sole exception of the deuteron. The D-wave decay widths of high spin, six quark states, such 
as the SIJP = 003+d* and the SIJp = -3^2+Nn, are in the range of tens of MeV to tens 
of keV and so these states might be narrow dibaryon resonances. 

The extended QDCSM is briefly introduced in Section II. In Section III, we present our 
results. The discussion and conclusion are given in Section IV. 



II. BRIEF DESCRIPTION OF THE EXTENDED QDCSM 



The QDCSM was put forward in the early 90's. Details can be found in Refs. J2|, |l2, 
|l3|. Although the short range repulsion and the intermediate range attraction of the A^A^ 
interaction are reproduced by the combination of quark delocalization and color screening, 
the effect of the long-range pion tail is missing in the QDCSM. Recently, the extended 
QDCSM was developed 7], which incorporates this long-range tail by adding vr-exchange 
but with a short-range cutoff to avoid double counting because the short and intermediate 
range interactions have been accounted for by the quark delocalization and color screening 
mechanism [ij]. The exchange of K and r] mesons has been shown to be negligible in this 
model approach y,|9|. Nevertheless, their effect, especially the tensor part, has been included 
in this calculation to check further whether they are negligible in our model approach. The 
extended QDCSM not only reproduces the properties of the deuteron well, but also improves 
agreement with NN scattering data as compared to previous work[l5|. 

The Hamiltonian of the extended QDCSM, wave functions and the necessary equations 
used in the current calculation are given below. The tensor interactions due to effective one 
gluon and the octet Goldstone boson exchanges are included. The details of the resonating- 
group method (RGM) have been presented in Refs. |lj, |l6[ . 

The Hamiltonian for the 3-quark system is the same as the well known quark potential 
model, the Isgur model. For the six-quark system, we assume 
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if i,j occur in the same baryon orbit, 
if i,j occur in different baryon orbits. 



otherwise, 

where tq is the short range cutoff for pion exchange between quarks. Qs is the vr quark 
couphng constant. ?7i,r is the measured vr mass. The K and rj meson exchange interactions, 
which have not been shown exphcitly in the above equation but have been included in this 
calculation, have a form very similar to that for the vr. The color screening constant, /i, is 
to be determined by fitting the deuteron mass in this model. All other symbols have their 
usual meanings, and the confinement potential Vcanfi^ij) has been discussed in Refs.|7l Il2|. 

The pion exchange interaction, V",r(^jj); affects only the u and d quarks. We take these 
to have a common mass, mg = rrid = rriu, i.e., ignoring isospin breaking effects. 

The quark wave function in a given nucleon (orbit) relative to a reference center (defined 
by S) is taken to have a Gaussian form characterized by a size parameter, b, 

The light quark mass, rrig, is chosen to be | of the nucleon mass. The strange quark mass, 
rris, baryon size parameter, b, effective quark-gluon coupling constant, a^, and the strength 
of confinement, Oc, are all determined by reproducing the nucleon mass, the A — N mass 
difference, an overall fit to other ground state baryon masses and by requiring the nucleon 
mass to be variational stable with respect to its size parameter, b. The quark-pion coupling 
constant g^ = ggq^, is obtained from the nucleon-pion coupling constant by a slight (< 10%) 
correction to the classic symmetry relation, viz., 

'^ = {M^lmgfQ''^e<'''\ (3) 

where Mat is the nucleon mass and the last factor provides the correction due to the extent 
of the quark wavefunction in the nucleon. The K and rj are assumed to have a flavor SU(3) 
symmetric quark-meson coupling constant and the same short range cutoff, tq, as the vr's. 
The color screening parameter, /i, has been determined by matching our calculation to the 
mass of the deuteron. All of the model parameters are listed in Table I. 



Table I: Model Parameters 










mq,ms{MeV) b{fm) a^{MeV ■ fm^"^) 


a. 


47r 


roifm) 


/i(/m-2) 


313,634 0.6022 25.03 


1.5547 


0.5926 


0.8 


0.90 



The model masses of all octet and decuplet baryons are listed in Table II. 



Table II: 


Single Baryon Masses in Units of MeV 










N S 


ASA 


S* 


'^* 


Q 


theor. 
expt. 


939.0 1217.5 
939 1193 


1116.9 1357.6 1232.0 
1116 1318 1232 


1359.6 

1385 


1499.7 
1533 


1652.3 
1672 



We use the RGM to carry out a dynamical calculation. The trial RGM di-baryon wave 
function is 

*(6g) =A[i^BA^,)^PB,U2)y'x{R), (4) 

where A is the antisymmetrization operator, ipsX^i) i = 1,2 is the baryon internal wave 
function including color-flavor-spin-orbital part, [■ ■ -Y^ means coupling the individual color- 
isospin-spin into the channel isospin-spin and overall color singlet. 

To simplify the RGM calculation, one usually introduces Gaussian functions with different 
reference centers Si i=l...n, which play the role of generating coordinates in this formalism, 
to expand the relative motion wave function x{R) of the two quark clusters, 

^^^) = ^iylf"T.C^e-'i^''-''^'''\ (5) 

In principle, any set of base wave functions can be used to expand the relative motion 
wave function. The choice of a Gaussian with the same size parameter, 6, as the single 
quark wave function given in Eq.((21), however, allows us to rewrite the resonating group 
wave function as a product of single quark wave functions; (see Eq.® below). This cluster 
wave function (physical basis) can be expressed in terms of the symmetry basis, classified 
by the symmetry properties in a group chain, which in turn allows the use of group th eory 
method to simplify the calculation of the matrix elements of the six quark Hamiltonian[l7[. 
In our calculations, we typically use 15 Gaussian functions to expand the relative motion 
wave function over the range 0-9 fm. For near threshold channels, such as the deuteron and 
H particle, more Gaussian functions are needed to extend the boundary to a larger extent 
to obtain more precise results as we have done previously. But in this calculation we did 
not make that effort because it is not necessary for our purpose. 

After including the wave function for the six-quark center-of-mass motion, the ansatz for 



the two-cluster wave function used in the RGM can be written as 

n „ 3 6 

j=l k Lfc=0,2 a=l /3=4 

where A; is the channel index. For example, for SIJ = —2, 0, 0, we have k = 1,2, 3, cor- 
responding to the channels AA, NS and SS. An angular momentum projection has been 
applied for the relative motion and L^ is the orbital angular momentum of the relative 
motion wave function of channel k. 

The delocalized orbital wavefunctions, ipnir. Si, e) and ipL^r, Si,e), are given by 

1 f Si Si 

i)R{r,Si,e) = ]^ l<^(^-y) + ^<^(^+y) 



N{e) = \/l + e2 + 2ee-^?/4^', 

where 0(^— y) ^^^ 0(^+ y) ^^^ ^^e single-particle Gaussian quark wave functions referred to 
above in Eq.(0), with different reference centers y and — y, respectively. The delocalization 
parameter, e, is a variational parameter determined by the djTiamics of the quark system 
rather than being treated as an adjustable parameter. The initial RGM equation is 

/ H{R, R')x{R')dR' = E f N{R, R')x{R')dR'. (8) 

With the above ansatz, the RGM Eq.® is converted into an algebraic eigenvalue equation, 

where N^j *' ' *", -f^j / '"' ' * are the wave function overlaps and Hamiltonian matrix elements, 
respectively, obtained for the wave functions of Eq.®. 

The partial width of a high spin dibaryon state decaying into a D-wave BB final state is 
calculated using Fermi's golden rule, in its nonrelativistic approximation, of course. Final 



state interactions have also been taken into account in our model approach |12i|. The decay 
width formula used in NVl -^ AS is a little different from that given in [l^l due to the fact 



that the decay products, A and S, are different particles with different masses, 



T{Nn ^ AE) = ^^ y -Atp V r WfAHn, (lo) 
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III. RESULTS 

Previously, we chose the di-fi as an example to study whether or not our model results 
were sensitive to the meson-exchange cut-off parameter, tq, and the result demonstrated that 
they are not[9|. Hence, we consider it is sufficient to calculate six-quark systems of different 
quantum numbers with a representative cutoff value of tq = O.Sfm. Table III displays the 
masses (in MeV) calculated for the dibaryon states of interest here. The deuteron channel 
result calculated previously is included in this table for comparison. It should be noted 
that in our calculation we assume the wavefunction to be zero at the boundary point, which 
is the usual boundary condition for bound states; i.e., we always solve the RGM Eq.Q 
as an eigenvalue problem. If the state is unbound, we will not obtain a stable minimum 
eigenenergy in the course of extending the boundary point. Therefore the "mass" listed 
in Table III for these unbound states is not a true mass of a dibaryon state. However the 
contributions of tensor interaction are still a meaningful measure. The masses listed under 
the MasSnt and Mass^t are the calculated masses of those channels without and with tensor 
interaction. 

Table HI : Tensor Interaction Effect on the Masses of Six-Quark Systems 
S,I,J channel Mass^i Mass^t 



0,0,1 NN 1880.295 1876.770(1880.250) 

-2,1,1 SS 2418.686 2418.517 

-3,1/2,1 AS 2477.491 2477.488 

-4,0,1 SS 2718.190 2718.189 

-5,1/2,1 EVL 3012.167 3012.140 



The first line of Table HI is the deuteron channel. If the tensor interaction is neglected the 
deuteron is unbound. Even if the tensor effect of gluon exchange is included, the deuteron 
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is still unbound. (The calculated mass is shown in parentheses). The tensor interaction due 
to vr exchange is critical to form the actual stable deuteron and the rj meson contribution 
is negligible. Fig. 1 shows the effective transition interactions of the NN S-D coupling due 
to gluon and vr exchanges. The rj contribution has not been shown because it is negligible. 
Obviously, the vr contribution is dominant. The boundary point is limited to 9 fm and 
coupling to AA channels has not been taken into account in this calculation; hence, the 
deuteron mass is a little higher than the best fit one which we reported previously. This 
weakness will not affect our conclusion regarding the tensor interaction effect mentioned 
above. 

The SS channel mass listed in the second line is lower than its own theoretical thresh- 
old 2435 MeV but higher than the NS, AS thresholds. A three S-wave channel coupling 
calculation has been done. The lowest mass is 2300 MeV which is still higher than the NS 
threshold and so can not form a narrow dibaryon resonance. The mass reduction due to the 
tensor interaction is very small as can be seen from Table III. Fig. 2 shows the S-D effective 
transition interactions in the SS channel due to vr, K and gluon exchange respectively. The 
IT contribution is reduced by about a factor of four. The gluon contribution is also very much 
reduced. In addition, it is repulsive and cancels the tt contribution. The K contribution is 
also repulsive and negligible. The rj contribution is even smaller. 

The AS, SS, and EQ states are near threshold. The tensor interactions in these channels 
are much less effective than that in the deuteron channel and not strong enough to bind these 
baryons into dibaryon resonances. The mass reductions due to the tensor interactions in 
these channels are negligible as shown in Table III. Fig. 3 shows the S-D effective transition 
interactions in the SS channel due to n, K and gluon respectively. The tt contribution is 
reduced even further than in the previous case. The gluon contribution is very much reduced 
in comparison with that in deuteron channel but a little bit enhanced in comparison with 
that in SS channel. The K contribution is further verified to be negligible and rj negligible 
also. 

Fig. 4 gives a direct comparison of the S-D effective transition interactions due to the 
tensor force of vr exchange in NN, SS and SS channels. Fig. 5 gives a direct comparison of 
the S-D effective transition interactions due to the tensor force of gluon exchange in NN, SS, 
AS, SS and Sf2 channels. The rj contribution has been calculated for all of these channels 
and all are negligible so they have not been shown explicitly. 

9 



As was mentioned in the introduction, the strangeness -3 LP = \2'^ Nfi has also been 
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predicted to be tightly bound j^ |10|, |11| if the tensor interaction is neglected. Taking into 
account the tensor force, this state can couple to D-wave AH and SH channels. Since AS is 
the lowest channel in all channels with strangeness S = —3, here we only take into account 
this channel. For the possible NQ{IJ = 1/2, 2) bound state we consider the two lowest AS 
D-wave decay channels: IS=l/2,0 and IS=1/2,1, respectively; here S specifies the channel 
spin. Such a tensor coupling has two effects: One is to modify the mass of the NQ state; 
the other is to induce a transition from the NQ to the D-wave AS final state and change the 
bound NQ to a resonance with finite width. Both of these effects have been calculated. 

The results show that, the mass of NQ{lJ=l/2,2) in the single channel approxima- 
tion is about 2566 MeV. Taking into account other S-wave channels coupling, such as 
SS, S*S, SS*, SAara(iS*S*, reduces the mass of the system to 2549 MeV, while adding the 
AS D-wave channel coupling only changes the value of mass slightly (not more than 1 MeV) . 

The widths of NQ decays to the AS D-wave with different spins are listed in Table IV. 
For comparison, the width of the d* decay to NN D-wave is also listed. 

Table IV. Decay width. I,S, J are isospin, spin and total angular momentum, respectively. 
Nn{IJ = 1/2, 2) -> AS D~ wave{S = 0,1 = 1/2, J = 2) F = U.OkeV 

Nn{IJ = 1/2, 2) -> AS D- wave{S =1,1 = 1/2, J = 2) F = 21.9keV 

d*{IJ = 0,3) ^NN D- wave{S = 1,1 = 0, J = 3) F = 6.57MeV 

From Table IV we see that the width of d* decay to the NN D-wave is 6.57 MeV. Com- 
parison with our previous results (l8l|. confirms that the width is not sensitive to the value of 
cutoff tq. The width of NQ decays to AS D-wave is about tens of keV, about three orders of 
magnitude smaller. The result is not changed significantly for decay channels with different 
spin. For example, the NQ -^ AS (spin=0) D-wave decay width is calculated to be 12.0 
keV, and the NQ —>■ AS (spin=l) D-wave decay width is 21.9 keV. These results confirm our 
expectation that the Nfi is a narrow dibaryon resonance. (The width of Nf2 -^ AS is smaller 
than that of the d* —>■ NN decay mainly due to the reduction of the tensor interaction, but 
also due to the fact that the N and A have the same flavor content while N, Q and A, S 
have differing flavor content in each baryon.) 
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IV. DISCUSSION AND CONCLUSION 

The effects of the tensor interactions of gluon and Goldstone boson exchanges on the 
dibaryon mass and decay width have been studied in the extended QDCSM. Only in the 
deuteron channel is the tensor interaction of vr exchange strong enough to bind the two 
nucleons into a loosely bound state. No other near threshold six-quark state studied so far 
in the u, d, s three flavor world can be bound together by the additional attraction induced 
by these tensor interactions. The S-D wave effective transition interactions due to vr and 
gluon tensor forces both decrease rapidly with increasing strangeness. In the SS and SS 
channels the effective transition interactions due to the gluon tensor term become repulsive 
and cancel the attractive vr contribution. The tensor contribution of explicit K and t] meson 
is confirmed as negligible due to the same short range truncation as for the tt, in our model 
approach y,|9| 

The mass shift of the IJ^ = ^2"'" NQ state induced by the tensor interaction is small (not 
more than 1 MeV) and the D-wave partial decay widths to AH with spin-0 and spin-1 are 
only 12.0 and 21.9 keV, respectively. Hence, the IJ^ = |2+ Nfi state appears to be a good 
candidate for a narrow dibaryon resonance. Altogether there are only two promising narrow 
dibaryon resonances in the light flavor world in our model approach: The IJ^ = 03^ d* and 



the IJP = i2+ Nfi. 



e in our model. However 
). Table II shows that 



The H particle and di-f2 are marginally strong interaction stab 
the theoretical binding energies of both are small (only few MeVJ4 
the calculated ground octet and decuplet baryon masses deviate from the measured ones 
about 18 MeV on average. A reasonable estimate of the model uncertainty for the dibaryon 
mass would be at least that large. Therefore, in our model, it is unjustified to assert 
that the H particle and di-fi might be strong interaction stable dibaryon candidates. This 
estimation is consistent with the latest di-A hypernuclear findings |l9l|. There are various 
broad resonances with widths ~ 150 -^ 250 MeV around the d* mass (~ 2180 MeV) in our 
model which makes the analysis of the NN scattering more difficult in the energy region 
2.1 —i> 2.4 MeV where a broad bump has been found in the pp and np total cross sections. 
We will report on those results later. The SIJ = —3, 1/2,2 Nf2 state is quite convincingly 
lower in mass than the NQ threshold, and quite possibly lower than the AHtt threshold, 
as well. We have shown the decay width to be as small as tens of keV. Such a narrow 
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dibaryon resonance might be detected by reconstructing the invariant mass of its two body 
decay products, A and S, in high Q production reactions using RHIC at Brookhaven and 
COMPASS at CERN and the future ones at JHF in Japan and FAIR in Germany. 

This model, the extended QDCSM, which proposes a new mechanism to describe the 
NN intermediate range attraction instead of the a meson, well describes, with the fewest 
parameters, the properties of the deuteron and the existing NN, NA and NS scattering data. 
Up to now, it is the only model which gives an explanation of the long standing fact that the 
nuclear and molecular forces are similar in character despite the obvious length and energy 
scale differences and that nuclei are well described as collection of A nucleons rather than 
3A quarks. In view of the fact that the H particle has not been observed experimentally, the 
BB interaction in the AA channel |l9l| predicted by this model may be a good approximation 
of the real world. Based on these facts we suppose the predictions about dibaryon states of 
this model might also be approximately correct. Of course, the QDCSM is only a model of 
QCD. The high spin, high strangeness dibaryon resonance, IJ^ = |2+ NQ, may be a good 
venue to search for new hadronic matter and to test whether or not the QDCSM mechanism 
for the intermediate range attraction is realistic. 

This work is supported by NSFC contracts 90103018 and 10375030. and by the U.S. 
Department of Energy under contract W-7405-ENG-36. F. Wang would like to thank the 
ITP for their support through the visiting program. 
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FIGURE CAPTIONS 

Fig.l The effective S-D wave transition interactions of gluon and n tensor force in the 
deuteron channel. 

Fig.2 The effective S-D wave transition interactions of gluon and vr, K tensor force in the 
IJP = 11+ SS channel. 

Fig. 3 The effective S-D wave transition interactions of gluon and vr, K tensor force in the 
Ijp = 01+ EE channel. 

Fig.4 A comparison of effective S-D wave transition interactions of vr tensor force in NN, 
SS, SS channels. 

Fig. 5 A comparison of effective S-D wave transition interactions of gluon tensor force in 
NN, SS, AS, SS, En channels. 
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FIG. 1: The effective S-D wave transition interactions of gluon and it tensor force in the deuteron 
channel. 
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FIG. 2: The effective S-D wave transition interactions of gluon and it, K tensor force in the 
IJP = 11+ SS channel. 
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FIG. 3: The effective S-D wave transition interactions of gluon and tt, K tensor force in the 
IJP = 01+ HH channel. 
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FIG. 4: A comparison of effective S-D wave transition interactions of vr tensor force in NN, SS, 
HH channels. 
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FIG. 5: A comparison of effective S-D wave transition interactions of gluon tensor force in NN, 

SS, AH, HH, Hi7 channels. 



19 



